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Abstract 

We show that the flux line lattice in a model layered high temperature su- 
perconductor becomes unstable above a critical magnetic field with respect 
to a plastic deformation via penetration of pairs of point-like disclination de- 
fects. The instability is characterized by the competition between the elastic 
and the pinning energies and is essentially assisted by softening of the lattice 
induced by a dimensional crossover of the fluctuations as field increases. We 
confirm through a computer simulation that this indeed may lead to a phase 
transition from crystalline order at low fields to a topologically disordered 
phase at higher fields. We propose that this mechanism provides a model of 
the low temperature field-driven disordering transition observed in neutron 
diffraction experiments on Eh^S^CaC^Os single crystals. 
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The mixed state in the high temperature superconductors(HTSC) has been studied from 
a wide variety of perspectives over the past few years |L|. Among the interesting issues under 
debates are the question of how vortices freeze in the so-called vortex glass phase |4| as 
well as the nature of the melted phase . In particular, one expects an intricate interplay 
of the underlying layered nature of the host material and the random point pins in shaping 
the nature of the frozen phase. 

A recent decoration performed on both sides of the Bi 2 Sr 2 CaCu208 single crystal sug- 
gests that vortex lines maintain their line integrity at least at very low field @. On the 
other hand, the rapid disappearence of Bragg peak intensity in the recent small angle neu- 
tron diffraction(SANS) [|J has been interpreted by the authors in terms of decomposition of 
lines into pancake vortices ||. A Lindemann-like criterion is usually invoked in discussing 
such a "decoupling theory" fLOf , but it fails to explain how the transition is brought about, 



if there is one. 

Traditionally, the role of a pinning potential has been treated as a perturbation within 
the harmonic elastic framework. In the weak pinning limit, the perturbed lattice tries 
to optimize its free energy by forming elastic domains of correlated region with minimum 
elastic energy while it is fragmented in larger length scales to take advantage of the random 



potential energy [|11|,|12| . Several attempts have been made recently to extend the simple 



dimensional argument |T3HT5| for length scales larger than the elastic volume, but without 
considerations of topological defects. Their effect is expected to dominate for strong enough 
disorder, resulting in a dislocation dominated " glass phase" |16|-p^]. 



In an earlier paper ||19|| , we showed that the dramatic features of the SANS results can not 
be accounted for by a dimensional crossover alone of a "clean" lattice. We further suggested 
that the flux lattice in the highly anisotropic B^S^CaC^Og may suffer an instability 
against penetration of topological defects and therefore may deteriorate rapidly across a 
characteristic flux line density. More recently, such a transition has been claimed to be 
observed also in a YBa 2 Cu307_,5 based superlattice ppfl . In this letter, we report computer 
simulation studies of this transition using simulated annealing on a model flux line system. 



We start by noting that when the density of flux line in a layered HTSC is varied, the effective 
anisotropy of the lattice and the disorder strength change. As a result, the elastic domains, 
described by an in-plane correlation length Rd and an out of plane correlation length Ld may 
shrink. In particular, when Ld reaches the interlayer spacing, a new low temperature glass 
state that is dominated by disclinations appears. We propose that the transition to this 
phase is first order and is characterized by an explosive invasion of point-like disclination 
pairs across a horizontal line in the B-T phase diagram. 



Employing a vortex representation of the Lawrence- Doniach model ||21|| , we consider a 
stack of coupled two dimesional vortex lattices with pancake vortex coordinates {r^} with 
i labelling the individual flux lines. An approximate pair-wise interaction under periodic 
boundary conditions has been derived to make large scale numerical calculations possible 
|22|| . With a choice of in-plane penetration depth A(0) = 1800A, d = lbX, k = 100 and 



anisotropy 7 = yM c /M q b = 55, we obtain, for a clean system, a melting line T m (B) which is 
in reasonable agreement with the experimental result |23| . The random potential is modelled 
by potential wells of uniform depth U p , of radius given by the lesser of 2£ ab and the grid size 
scattered at random positions of each layer with an areal density of n p = l/a p = B p /(j) . 
We choose U p = 5 87T 2x2^ , B p = 150G in 16(Set #1) and 32(Set #3) layers which give the 
freezing temperature ~ 40°K for lkG < B < 20kG under a Lindemann-like criterion in 
reasonable agreement with the experimentally determined irreversibility line T irr (B). We 
also have results for a different pin density B p = 2AkG in 16 layers. (Set #2) The potential 
used in our simulation may be viewed as a coarse grained potential fluctuations with an 
unspecified microscopic origin. The line density(100G ~ 2kG) was varied by changing sizes 
of grid cells and pin densities for a fixed number of vortices(64 lines in L = 16 or 32 layers) 
22| . To obtain low temperature properties, we use a simulated annealing procedure starting 
with a perfect triangular lattice at T = 1.3 ■ Tf(B) ~ 50°-fT and gradually decreasing T 
with steps of dT = 5°K over 40,000 Monte Carlo steps. At T = 4°K, additional 20,000 
steps were performed to make measurement of physical quantities. This formally resembles 
a field-cooling procedure employed in typical measurements. The number of different pin 



configurations was limited to five for each line density for practical reasons. The disclination 
charge density n^ii, z) is measured every 50 steps through Delaunay triangulation performed 



in each layer to determine the coordination number Z of each vortex [23]. Then a charge of 
q = (Z — 6) is assigned to rid(i, z). 

Results In the inset to Fig. 1, we show the partial Fourier transform S(q x , q y , z = L) of 
the vortex density correlation function / d 2 p exp[iq-p] < n v (p, z = L)-n v (0, 0) > as evaluated 
in our simulation where n v (p, z) is the local vortex density. Ideally, SANS with q z ~ aims 
to measure / dzS(q x ,q y , z) and therefore for a macroscopic sample(L — > oo), we expect it 
to be related to the diffraction pattern of ||. The simulated Bragg peak intensities(dots 
in Fig. 1), decaying rapidly over a narrow range of 5B/B CZ ~ 0.3, looks strikingly similar 
to the appropriately scaled data from Cubitt et al. The values of B cz were 1800(Set #1) 
and 550(^2) Gauss while we used 500 Gauss for the experimental data. The statistical 
deviation of the intensity, SIgi/Igi, over five disorder configurations is within 20 %. The 
nature of the apparent peak for Set #1 in the figure, therefore, can not be resolved at this 
time. To verify what is driving the rapid drop, we look into the behavior of topological 
defects and the fraction of pinned vortices, f p which is given by (j2i, z @( r i,z)) /N where 
0(r) gives 1 if r is within any of the pinning wells and otherwise. In Fig. 2, we observe a 
jump in f p within the narrow range of SB/B CZ ~ 0.2, the same region over which the Bragg 
intensity vanishes. There is a corresponding increase in the total number of disclinations, 
J2i,z i\ n d(h z )\) accompanied by a rapid drop of the in-plane hexatic order parameter, \1> 6 
|25|| . As shown in the inset to the low field side, a typical configuration has a finite number 
of defects, but they in general appear as neutral disclination pairs or quartets(equivalent 
to a bound pair of dislocations), healing each other over a finite length along z-axis. With 
these bound defects, the lattice order is disrupted only over a finite distance given by the 
size of these bound defects, and manifest itself as a distinct Bragg spots of Fig. 1. It is to 
be contrasted with the higher field configuration in which defects are threading from-top-to- 
bottom layers of the sample in various sizes and separations. In this phase, we find that the 
weak hexatic order is still present in each layer and that their correlation along the c-axis is 



sensitively dependent on the overall strength of the pins |26[ . 

Introducing an integer variable n d (i, z) which gives 1 for a disclination and otherwise 
for lattice site (i,z), we define the probability distribution for length I of continuous defect 
lines: 

z +l 

nO=[(££[l-<(Mo)][l-r4(^o + m)] II n* d (i,z))], (1) 

i,Z Q 2 = 20 + 1 

where < . . . > and [. . .] mean thermal and defect averages, respectively. In a pin-free system, 
the distribution above T m (B) shows a gradual crossover behavior as B increases, reflecting 



the underlying dimensional crossover of the fluctuations of the lattice |23| , |27l1 . In this case 
of T <C Tf with quenched disorder, however, we observe a more dramatic change which 
suggests a field-driven phase transition. In Fig. 3, we show V(l) of continuous defect length 
I for fields B/B cz = 0.2 ~ 2 averaged over five different pin configurations (Set #1). Despite 
large statistical fluctuations at the large tail, we can clearly see that V{1) for high field ~ e~ al 
qualitatively differs from the low field distribution for which V(l) = for I > l max (B). The 
exponential dependence for high field can easily be understood in terms of defects generated 
independently in each layer with a probability p. A probability for an accidentally aligned 
line of length I will then be V(l) ~ (1 —p) 2 p l / d ~ exp[— ^-| logp|]. In the low field, the length 
of defect line is truncated by the energetic balance between the defect rigidity and the overall 
gain from random pins. 

This observation gives a clue in understanding the role of pins in the defect representation 
of the free energy. Clearly, the pinning potential acts as an effective "temperature" to drive 
the topological defects into the lattice by encouraging vortices to make large excursions to 
seek optimal pinning configurations. The occurrence of these topological defects signals the 
breakdown of the elastic theory and the transition to a new defect-dominated phase. We 
can estimate the breakdown field B cz as follows: Let us assume that the accumulated vortex 
displacements Ui(x) reach 0(as) over a volume Vd of radius Rd in the afe-plane, length Ld 
along the c-axis. For B < [Up/ed^B 2 / Bj, where Bj = (po/^d) 2 , a variational minimization 
of elastic free energy yields Rd/ciB < 1, suggesting a breakdown of the quasi- crystalline order 



occurring in the low field region of the phase diagram. Consequently, we assume a very short 
ranged in-plane order and fix ~ as, and perform a partial variational calculation for a 
tube of variable length Ld\R d ~ aB to obtain, Ld(B)/d ~ [e 2 d 2 B 2 j/U 2 B p B} 1 ^ . The breakdown 
field B cz is obtained from the condition Ld(B cz ) < d. Here we assume that the statistical gain 
in pinning energy by bending of such tubes is given by U p (n p Vd/ d)°' 5+s / (a 2 B L d ) with 5 = 0, 
a rough approximation. Using the parameter values we chose, this yields B cz ~ 2.4kG,[Set 
#1] in reasonable agreement with the result of the simulation. 

The results of the simulations may be summarized in a phase diagram given in Fig. 4. 
At low temperatures, the effective "noise temperature" of the random pinning field drives 
the system from a quasi-crystalline phase at low fields to a defect-dominated topologically 
disordered phase at higher fields. This phase transition is facilitated by weak interlayer 
coupling, thus explaining the small value of the critical field B cz for highly anisotropic 
materials such as E^S^CaC^Os . On the other hand, for less anisotropic materials such as 
YBa 2 Cu 3 07_ ( 5 , this instability should be pushed to a higher field. Because the disorder is 
induced by the pins, the vortices are still highly localized, hence in a glassy phase. We have 
not yet used our simulations to study the kinetics of this phase, however, from our previous 
studies we expect the kinetics of the phase to be sub-ohmic, so that it should still be 
a superconductor. 

As the temperature is increased in the topologically disordered phase, we expect to 
reach a "depinning temperature" above which the defects become liquid-like and the trans- 
port would be ohmic. The exact shape of the line of phase transitions between the quasi- 
crystalline phase and the topological glass phase depend on the details of the competing 
energies and may terminate at the melting line T^f(B) which was recently shown to bend 
toward and terminate at a field of ~ 500 Gauss |[29| . In the vicinity of the transition line 
B cz , the explosion of topological defects will be expected to dominate the kinetics and may 
provide an explanation of the "fishtail" peak-effect in the critical current observed as the 
B-field is varied in this region of the phase diagram P0|PT[ ]. The results reported here are 
reminiscent of a transition observed independently by Gingras and Huse from simulations 
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of a ferromagnetic 3d XY model with a random field [32 
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FIGURES 

FIG. 1. The relative integrated intensity of the first order Bragg peak in the partial structure 
factor for Set #1 and #2. The intensity was normalized by that of the clean sample at the same 
temperature and field. The broken line is from the Neutron scattering data Q drawn with an 
arbitrary scale. The field was arbitrarily scaled by a value at which the intensity reaches the flat 
bottom. The insets show the simulated diffraction pattern for three representative points (labelled 
a,b,c) of Set #1 in the curve. Intensity for (c) was multiplied by a factor of 30. 

FIG. 2. Pinned fraction of vortices across the instability line for Set #1. The broken line is a 
guide showing the naive trend of ~ B p /B as expected from the relative ratio of the densities. The 
abrupt increase in the pinned fraction is closely accompanied by the proliferation of topological 
defects as shown as black(z > 6) and gray(z < 6) dots in the insets. 

FIG. 3. V(l) for various values of B averaged over five different realizations of pins(Set #3). 
Number of layers was 32. Filled circles correspond to the hexatic glass region while open symbols 
are for higher field topological glass regime. Broken lines are guides for the eye. 

FIG. 4. Schematic phase diagram with strong point pins. The low field elastic lattice is turned 
into a topological glass state as B increases over B cz . The freezing temperatures (Tf d and Tj d ) are 
characterized by the dimensionality of topological defects as discussed in the text. We also note 
that there is a distinct decoupling crossover line Td c in the liquid phase across which the vortex 
lines disintegrate into pancake vortices through cutting and reconnection [ 33 1 . 
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